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Aaronson and Arkhipov recently used computational complexity theory to argue that classical 
computers very likely cannot efficiently simulate linear, multimode, quantum-optical interferometers 
with arbitrary Fock-state inputs [S. Aaronson and A. Arkhipov, Theory of Computing, 9, 143 
(2013)]. Here we present an elementary argument that utilizes only techniques from quantum 
optics. We explicitly construct the Hilbert space for such an interferometer and show that that 
its dimension scales exponentially with all the physical resources. Finally we also show in a simple 
example just how the Schrodinger and Heisenberg pictures of quantum theory, while mathematically 
equivalent, are not in general computationally equivalent. 



There is a history of attempts to use linear quan- 
tum interferometers to design a quantum computer. 
Cerny showed that a linear interferometer could solve 
NP~complete problems in polynomial time but only with 
an exponential overhead in energy [jj] . Clauser and Bowl- 
ing showed that a linear interferometer could factor large 
numbers in polynomial time but only with exponential 
overhead in both energy and spatial dimension [2| . Cerf, 
Adami, and Kwiat showed how to build a programmable 
linear quantum optical computer but with an exponential 
overhead in spatial dimension [3| . 

Nonlinear optics provides a well-known route to uni- 
versal quantum computing |4<] . We include in this nonlin- 
ear class the so-called 'linear' optical approach to quan- 
tum computing [5|, because this scheme contains an ef- 
fective Kerr nonlinearity [6J. 

In light of these results there arose a widely held belief 
that linear interferometers alone, even with nonclassical 
input states, cannot provide a road to universal quan- 
tum computation and, as a corollary, that all such de- 
vices can be efficiently simulated classically. However, 
recently Aaronson and Arkhipov (AA) gave an argument 
that multimode, linear, quantum optical interferometers 
with arbitrary Fock-state photon inputs very likely could 
not be simulated efficiently with a classical computer 7] . 
Their argument, couched in the language of quantum 
computer complexity class theory, is not easy to follow 
for those not skilled in that art. Nevertheless, White, 
collaborators, and several other groups carried out ex- 
periments that demonstrated that the conclusion of AA 
holds up for small photon numbers @t11| ■ Our goal here 
is to understand — from a physical point of view — why 
such a device cannot be simulated classically. 

We independently came to the same conclusion as AA 
in our recent analysis of multi-photon quantum random 
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FIG. 1. The quantum pachinko machine for numerical depth 
1/ = 3. We indicate an arbitrary bosonic dual-Fock in- 
put |A'^) I A/) at the top of the interferometer and then the 
lattice of beam splitters (B), phase shifters (</?), and pho- 
ton-number-resolving detectors (V). The vacuum input 
modes |0) (dashed lines) and internal modes |V') (solid lines) 
are also shown. The notation is such that the superscripts la- 
bel the level i and the subscripts label the row element from 
left to right. 



walks in a particular multimode interferometer called a 
quantum "pachinko" machine shown in FiglT] [I2l]. The 
dual-photon Fock state |A'') \M) is inputted into the top 
of the interferometer and then the photons are allowed 
to cascade downwards through the lattice of beam split- 
ters (B) and phase shifters {(p) to arrive at an array of 
photon-number-resolving detectors (V) at the bottom. 
Our goal (that of a pachinko player) was to calculate 
the probability that p photons were detected at the q*^ 
detector for arbitrary input photon number and lattice 
dimension. We failed utterly. It is easy to see why. 

Working in the Schrodinger picture, we set out to com- 
pute the probability amplitudes at each detector by fol- 
lowing the Feynman-like paths of each photon through 
the machine, and then summing their contributions at 
the end. For a machine of numerical depth L, as shown 



in Fig. [1] it is easy to compute that the number of such 
Feynman-hke paths is 2^(^+*^). So for even a meager 
number of photons and levels the solution to the problem 
by this Schrodinger picture approach becomes rapidly in- 
tractable. For example, choosing N = M — 9 and L — 6, 
we have 2^** = 5 x 10^^ total possible paths, which is 
about four orders of magnitude larger then the number 
of atoms in the observable universe. We were puzzled by 
this conclusion; we expected any passive linear quantum 
optical interferometer to be efficiently simulatable classi- 
cally. With the A A result now in hand, we set out here 
to investigate the issue of the complexity of our quantum 
pachinko machine from an intuitive physical perspective. 
The most mathematics and physics we shall need is ele- 
mentary combinatorics and quantum optics. 

Following Feynman, we shall explicitly construct the 
pachinko machine's Hilbert state space for an arbitrary 
level L, and for arbitrary photon input number, and show 
that the space's dimension grows exponentially as a func- 
tion of each of the physical resources needed to build 
and run the interferometer [13| . From this result we 
conclude that it is very likely that any classical com- 
puter that tries to simulate the operation of the quan- 
tum pachinko machine will always suffer an exponential 
slowdown. We will also show that no exponential growth 
occurs if Fock states are replaced with photonic coherent 
states or squeezed states, which elucidates the special 
nature of photonic Fock states. 

As our argument is all about counting resources, we 
have carefully labeled all the components in the pachinko 
machine in Fig. [T]to help us with that reckoning. The 
machine has a total of L levels of physical depth d 
each. The input state at the top is the dual-Fock state 
\N)^ IM)^, where the superscripts label the level num- 
ber and the subscripts the element in the row at that 
level (from left to right). We illustrate a machine of to- 
tal numerical depth of L = 3. For 1 < ^ < L, we show 
the vacuum input modes along the edges of the machine. 
The resources we are most concerned about are energy, 
time, spatial dimension, and number of physical elements 
needed to construct the device. All of these scale either 
linearly or quadratically in either L or A^ -1- Af . The to- 
tal physical depth is D = Ld and so the spatial area 
is A = {'\/2DY = 2L'^d^. Using identical photons of 
frequency w, the energy per run is E — [N + AI)huj. 
The time it takes for the photons to arrive at the de- 
tectors is T = ^/2Ld/c, where c is the speed of light. 
In each level the photons encounter i number of BS so 
the total number is #B = J2i=i ^ = L{L -t- l)/2. Be- 
low each BS (with the exception of the L*'' level) there 
are two independently tunable PS for a total number of 
PS that is #(^ = Xlfci^ 2£ = L{L - 1). The total num- 
ber of detectors is ^^V = 2L. The total number of in- 
put modes is equal to the total number of output modes 
and is #/ = #0 = 2L. The total number of internal 
modes is #V = Sf=i ^^ = L{L — 1). As promised ev- 
erything scales either linearly or quadratically in either 
L or N + M. 



The input state may be written in the Heisenberg 
picture as \N)\\M)l = (aj°)^(4°)*^ |0)" 10)° /VaHM!, 
where a^ is a modal creation operator. Each BS per- 
forms a forward unitary mode transformation, which we 
illustrate with B\, of the form d} = irjdj -I- tja^ and 
a\ = t\a\ + ir\d^ where the reflection and transmission 
coefficients r and t are positive real numbers such that 
r"^ Jr t^ = R + T ^ I. The choice r = t = 1/^2 im- 
plements a 50-50 BS. Each PS is implemented by, for 
example, applying the unitary operation exp(i(p}n}) on 
mode |'0)i; where n\ :— d\ d\ is the number operator, d\ 
is the annihilation operator conjugate to a\ , and (p\ is 
a real number. Finally the 2L detectors in the final level 
L are each photon number resolving 14]. 

To argue that this machine (or any like it) cannot 
be simulated classically, in general, it suffices to show 
that this is so for a particular simplified example. We 
now take N and L arbitrary but M — and turn off 
all the phase shifts and make all the BS identical by 



setting ifil = 0, tl 



t, and 



r for all (kj). 



We then need the backwards BS transformation on the 



creation operators, which is, 
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i| = ira[ + td2 and 
Similar transforms apply down the 
machine at each level. With Af — the input simpli- 
fies to \N)° |G)° ^ (dj°)^ |G)° |0)° /VM and now we ap- 
ply the first backwards BS transformation \ijj)i \'4')2 — 
{ira[ 
one 
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+ tal^)'^ |0)i 10)2 /\/M to get the state at level 



At every new level each d^ will again bifurcate accord- 
ing to the BS transformations for that level, with the 
total number of bifurcations equal to the total number 
of BS, and so the computation of all the terms at the 
final level involves a polynomial number of steps in L. It 
is instructive to carry this process out explicitly to level 
i = 3 to get, 
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where we have used a tensor product notation for the 
states. If r = or r = 1 the state is easily computed. 
Since we are seeking a regime that cannot be simulated 
classically we work with r = t = l/\/2. 

It is now clear from Eq.([T|) what the general form of 
the solution will be. We define 
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and the general solution has the form, 
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combinatorics and is: 



\^)'-^,(t^f^') 



\or 



m 



E 






N 
ni,n2,. ■ ■,n2L 

L 



(3) 



l<k<2L 



where all the coefficients af will be nonzero in general. 
Since all the operators commute, as they each operate 
on a different mode, we have expanded Eq.Q using the 
multinomial theorem where the sum in the expansion 
is over all combinations of non-negative integers con- 
strained by A'^ = '^£=1 '^i ^^'^ 
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is the multinomial coefficient [lal- The state \tp) is 
highly entangled over the number-path degrees of free- 
dom. Each monomial in the expansion of Eq. ^ is unique 
and so the action of the set of all monomial opera- 
tors on the vacuum will produce a complete orthonor- 
mal basis set for the Hilbert space at level L, given by 
IV') f„ } •= 11^=1 l"-^)f J where the ni are subject to the 
same sum constraint. Let us call the dimension of that 
Hilbert space dim.[H{N,L)], which is therefore the total 
number of such basis vectors. 

Taking L = 3 and TV = 2, we can use Eq.® to com- 
pute the probability a particular sequence of detectors 
will fire with particular photon numbers. What is the 
probability detector one gets one photon, detector two 
also gets one, and all the rest get zero? This is the mod- 
ulus squared of the probability amplitude of the state 
|1>? |1>2 10)3 10)4 10)5 10)6- Setting r = t= 1/^2 for the 
50-50 BS case, from Eq.dl]) we read off af — irt^ — 
i/{2y/2) and al = -r'^t = -l/(2%/2), and so the proba- 
bility of this event is given by Pnoooo — 0.031. 

It turns out that it is possible (for general L and N) 
to compute the binary joint probabilities, that detector 
p gets n photons and detector q gets m [I^. However 
the pachinko player instead wants to know which de- 
tector gets the most photons on average. For L = 3 
and A^ = 2, to compute this he must calculate all the 
dim[iJ(2, 3)] = 21 probabilities Pnpqrst ^'^d then choose 
from that list which detector's likely to detect the most 
photons. However for arbitrary L and A^ this optimiza- 
tion he can not do, as we now show. We can provide 
a closed form expression for dim[_ff (iV, L)] by realizing 
that it is the same as the number of different ways one 
can add up non-negative integers that total to fixed N . 
More physically this is the number of possible ways that 
N indistinguishable photons may be distributed over 2L 
detectors. The answer is well known in the theory of 



dvm\H{N, L)] = 



N + 2L-1 

N 



(5) 



where this is the ordinary binomial coefficient |17| . For 
our example with L — 3, N = 2, Eq.® implies that the 
number of distinct probabilities Pnp„rst to be tabulated 
is again 21. 

We first examine two "computationally simple" ex- 
amples. Taking TV arbitrary and L — 1 we get 
diin[H{N, 1)] — TV -I- 1, which is easily seen to be the 
number of ways to distribute N photons over two de- 
tectors. Next taking A'^ = 1 and L arbitrary we get 
dim[iJ(l, L)] — 2L, which is the number of ways to dis- 
tribute a single photon over 2L detectors. If we were 
to invoke Dirac's edict — "Each photon then interferes 
only with itself." — we would then expect that adding 
a second photon should only double this latter result 
[18l |. Instead the effect of two-photon interference on 
the state space can be seen immediately by computing 
dim[iJ(2, L)] = L{2L + 1). That is adding a second pho- 
ton causes a quadratic (as opposed to linear) jump in the 
size of the Hilbert space. Dirac was wrong, photons do 
interfere with each other, and that multiphoton interfer- 
ence directly affects the computational complexity. All 
these three cases are simulatable in polynomial time steps 
with A'^ and L, but we see a quadratic jump in dimension 
as soon as we go from one to two photons. These jumps 
in complexity continue for each additional photon added 
and the dimension grows rapidly. 

We therefore next investigate a "computationally com- 
plex" intermediate regime by fixing A^ — 2L — 1. 
That is we build a machine with total number of lev- 
els L and then choose an odd-numbered photon in- 
put so that this restriction holds. Eq.© becomes 
dim[iJ(A^)] = {2Ny./{N[f. Deploying Sterling's approx- 
imation for large A^, in the form n! = (n/e)"\/27rn we 
have dim[i?(A^)] = 2'^^ / VnN . This is our primary re- 
sult. The Hilbert space dimension scales exponentially 
with A^ = 2L — 1. Since all the physical parameters 
needed to construct and run our quantum pachinko ma- 
chine scale only linear or quadratically with respect to A^ 
or L, we have an exponentially large Hilbert space pro- 
duced from a polynomial number of physical resources 
— Feynman's necessary condition for a potential univer- 
sal quantum computer. 

Let us suppose we build onto an integrated optical cir- 
cuit a machine of depth L = 69 and fix A^ = 2L—1 = 137. 
Such a device is not too far off on the current quantum 
optical technological growth curve [l^- Then we have 
dim[iJ(137)] == 10^^, which is again on the order of the 
number of atoms in the observable universe. Following 
Feynman's lead, we conclude that, due to this exponen- 
tially large Hilbert space, it is likely that a classical com- 
puter can not efficiently simulate our quantum pachinko 
machine, in particular, and hence it likely can not in 
general efficiently simulate all multimode linear optical 
interferometers with arbitrary Fock-state inputs. 



Let us now compare our Heisenberg picture result to 
that of the Schrodinger picture. In the computationaUy 
complex regime where N = 2L — 1 the number of distinct 
Feynman-like paths we must follow in the Schrodinger 
picture is 2^^ = 2^(^+1)/^ ^ 2^'^\ Taking N = 137 
and L = 69, as in the previous example, we get an as- 
tounding 2^"'^'^ = 4 X 10^^^^ total paths. Dirac proved 
that the Heisenberg and Schrodinger pictures are math- 
ematically equivalent, that they always give the same 
predictions, but we see here that they are not always 
necessarily computationally equivalent [20|. Calculations 
in the Heisenberg picture are often much simpler than 
in the Schrodinger picture. The fact that the two pic- 
tures are not computationally equivalent is implicit in 
the Gottesman-Knill theorem, however it is satisfying to 
see here just how that is so in a simple optical interfer- 
ometer |2l| . 

To contrast this exponential overhead from the re- 
source of bosonic Fock states, let us now carry out the 
same analysis with the bosonic coherent input state in- 
put 1/3)]^ 10)2, where we take the mean number of photons 
to be |/3p = n. In the Heisenberg picture this input be- 
comes D?(/3) \0)1 10)2 where_^(/3) = exp(/3al°-/3*a?) is 
the displacement operator [22| . Applying the BS trans- 
formations down to final level L we get. 
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At the output we have 2L coherent states that have been 
modified in phase and amplitude. This is to be expected, 
as it is well known that linear interferometers transform 
a coherent state into another coherent state. Since all 
the coefficients af are computable in #5 ~ L{L + l)/2 
steps, this result is obtained in polynomial time steps in 
L, independent of n. The mean number of photons at 
each detector is then simply nf = \j3a^\'^ = n\aj'\'^. 

A similar analysis may be carried out for bosonic 
squeezed input states. Taking, for example, a 
single-mode squeezed vacuum input |C)i |0)2 = 
Si{^) |0)^ 10)2, with the squeezing operator defined as 
5?(?) = exp[(r(a;)' - a&\y)m, we arrive at. 
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which does not in general decompose into a separable 
product of single-mode squeezers on each output port. 
Nevertheless the probability amplitudes may still be com- 
puted in a time polynomial in L by noting that, from 
Eq.(l5]) with N = 2, there are at most 2L{L + 1) terms 



in this exponent that must be evaluated. This result 
generalizes to arbitrary squeezed-state inputs Q. The 
output of the interferometer may be then calculated on 
the transformed device in polynomial steps in L. 

The exponential scaling comes from the bosonic Fock 
structure |A^) = (a^)^ |0) /y/W. and the fact that N can 
be made close to L with identical bosons in a single 
mode. It is well known that beam splitters can gener- 
ate number-path entanglement from separable bosonic 
Fock states. For example, the simplest version of the 
Hong-Ou-Mandel effect at level one with separable in- 
put |l)j 11)2 becomes \ip)^ \ip)2 = (ia} -I- al ){a\ + 
iajl) \0)l 10)2 /2 = i[\2)\ |0)2-l- 10)} |2)2]/%/2 a NOON state 
[24| . Such entangled NOON states violate a Bell inequal- 
ity and are hence nonlocal even though the input was not 
[25| . For arbitrary bosonic Fock input states and inter- 
ferometer size the amount of number-path entanglement 
grows exponentially fast. 

In conclusion, we have shown that a multimode linear 
optical interferometer with arbitrary Fock input states 
is very likely not simulatable classically. Our result is 
consistent with the argument of AA. Without invoking 
complexity theory, we have argued this by explicitly con- 
structing the Hilbert state space of a particular such in- 
terferometer and showed that the dimension grows expo- 
nentially with the size of the machine. The output state 
is highly entangled in the photon number and path de- 
grees of freedom. We have also shown that simulating the 
device has radically different computational overheads in 
the Heisenberg versus the Schrodinger picture, illustrat- 
ing just how the the two pictures are not in general com- 
putationally equivalent within this simple linear optical 
example. It is unknown if such bosonic multi-mode in- 
terferometers as these are universal quantum computers, 
but regardless they will certainly not be fault tolerant. 
As pointed out by Rohde [26| , it is well known that Fock 
states of high photon number are particularly sensitive 
to loss [2711 . They are also super-sensitive to dephasing 
as well [28|. This implies that even if such a machine 
turns out to be universal it would require some type of 
error correction to run fault tolerantly. Nevertheless such 
devices could be interesting tools for studying the rela- 
tionship between multiphoton interference and quantum 
information processing for small numbers of photons. If 
we choose each of the PS and BS transformations inde- 
pendently of each other, we have a mechanism to program 
the pachinko machine by steering the output into any of 
the possible output states. Even if universality turns out 
to be lacking we may very well be able to exploit this pro- 
grammability to make a special purpose quantum simu- 
lator for certain physics problems such as frustrated spin 
systems [29{ . 
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